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We use a quantum trajectory-based semi-classical method to account for Coulomb interaction 
between the photoelectron and the parent ion in the classically forbidden, sub-barrier region during 
strong-field tunneling ionization processes. We show that — besides the well-known modification 
of the tunneling ionization probability — there is also an influence on the interference pattern in 
the photoelectron spectra. In the long- wavelength limit, the shift of the intra-cycle interference 
fringes caused by sub-barrier Coulomb effects in the laser polarization direction can be derived 
analytically. We compare our results with ab initio solutions of the time-dependent Schrodinger 
equation and find good agreement in the long-wavelength regime, whereas the standard strong field 
approximation fails. We show that the nodal structure along low-order above-threshold ionization 
rings is also affected by sub-barrier Coulomb effects. 
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I. INTRODUCTION 

It has been shown experimentally that photoelectron 
spectra obtained from strong-laser matter interaction 
contain a wealth of information about the driving laser 
field and the target [H-Ql. Decoding this information and 
discriminating the various spectral features with respect 
to laser parameters, target geometry, many-electron ef- 
fects etc. is often difficult. On one hand, experiments 
and ab initio solutions of the time-dependent Schrodinger 
equation (TDSE) yield the "exact" spectra but little in- 
sight. On the other hand, the widely used workhorse of 
strong-field physics, the strong field approximation (SFA) 
, allows for an intuitive interpretation when formulated 
in terms of interfering quantum orbits via the saddle- 
point approximation (see, e.g., However, the SFA 
is rather inaccurate for the case of long-range Coulomb 
potentials due to the neglect of the interaction between 
the photoelectron and the parent ion after the time of 
ionization. In fact, various spectral features originating 
from this Coulomb interaction have been discovered re- 
cently among them the so-called "low-energy structure" 
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asymmetries in elliptically polarized laser fields 
13l | , and holographic side-lobes [1, 01 ■ 
Within the trajectory-based SFA, the matrix element 
for ionization is written as a coherent sum over ioniza- 
tion times for which the corresponding quantum orbits 
all lead to a given, final momentum p at the detector 
(see, e.g., If there are more than one such or- 

bits of similar weight, interference comes into play. In- 
terfering orbits originating from subsequent laser cycles 
give rise to peaks in energy spectra that are separated 
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by TiLO with uj being the laser frequency, the so-called 
"above-threshold ionization" (ATI) peaks. We call this 
"inter-cycle interference", in contrast to "intra-cycle in- 
terference" (l^ . [l5| where the interfering orbits originate 
from within one laser cycle, leading to, e.g., the holo- 
graphic interference patterns observed in [H, Intra- 
and inter-cycle interference in time is in analogy to a 
double slit and a grating in position space, respectively. 

Once the SFA is formulated in terms of trajectories, 
one may Coulomb-correct them 0, [l^ Hi [Hj [13 ■ This is 
straight-forwardly done in the classically allowed region, 
i.e., from the so-called "tunnel exit" to the detector. In 
the current paper we will extend the Coulomb-correction 
of photoelectron spectra to the sub-barrier region where 
time and orbit are complex. The imaginary time method 
has since many years been employed for the derivation 
of Coulomb-corrected tunneling ionization rates (see (l8| 
for a review). In fact, the difference in the probability for 
tunneling through a triangular barrier of a short-range 
potential and through the corresponding Coulomb bar- 
rier can be many orders of magnitude (l9l. [20j. However, 
to the best of our knowledge sub-barrier Coulomb effects 
on the interference pattern of photoelectron spectra have 
not been investigated. 

The outline of our paper is as follows. In Section HIl we 
review the trajectory-based Coulomb-corrected strong- 
field approximation (TCSFA) in real time and space be- 
fore we introduce, in Subsection III Bl the sub-barrier 
Coulomb correction (sub-CC). In Section Hill results on 
the intra-cycle interference are presented. In particular, 
the fringe shift in polarization direction due to the sub- 
CC is derived and discussed. There is no sub-CC ex- 
actly in perpendicular direction, as discussed in Subsec- 
tion |IIIB] The nodal structure along low-order ATI rings 
is investigated in Section lTVl We summarize in SectionlVl 
Atomic units h = rUf, = lei = 47reo = 1 are used, unless 
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noted otherwise. 



II. TRAJECTORY-BASED 
COULOMB-CORRECTED STRONG-FIELD 
APPROXIMATION (TCSFA) 

A quantum orbit in plain SFA is a semi-classical tra- 
jectory of a photoelectron starting at a complex time tg , 
subsequently governed by the external light field, and fi- 
nally detected with the asymptotic momentum p. The 
transition amplitude between the bound state to the con- 
tinuum state of asymptotic momentum p is a coherent 
superposition of contributions from all relevant quantum 
orbits (indexed by a) that lead to the asymptotic mo- 
mentum p, 
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Cp{ti"^) is a prefactor whose form is not of interest in 
this work (see, e.g., 0, IMI for more details). 



is the action in which 

v{t) =p + A{t) 



dt 



(2) 



(3) 



with the vector potential of the laser field A{t) in dipole 
approximation, and Ip is the ionization potential. All the 
trajectories a leading to an asymptotic electron momen- 
tum p fulfill the equation of motion 



dv 
It 



= -E{t), E{t) = -dtA{t) 



(4) 



and start at the complex saddle-point (i.e., ionization) 

time ti"' , which is determined by the saddle point equa- 
tion (SPE) 

d 



dt 



Wo{t) 



= 



p + A{t 



With Ip > and p real-valued, ti"-* is complex. From 
(j4|) and ([3]) follows that the canonical momentum is con- 
served, i.e., p = const. The amplitude ([T]) clearly demon- 
strates the multi-slit in time-nature of quantum mechan- 
ical ionization dynamics Q and lays the cornerstone for 
our further analysis. 

The key idea of the TCSFA is to insert the Coulomb 
force due to the ion of nuclear charge Z into the electronic 
equation of motion. 
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dv „, , Zr 
— = -E(t) - —. 
dt ^' rP' 



(6) 



and the corresponding Coulomb potential into the action. 
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Figure 1. (color online), (a) Integration path in the complex 
time plane. Starting at ts we integrate parallel to the imag- 
inary time axis down to tr on the real axis, and then along 
the real axis to infinity. The corresponding situations in po- 
sition space around the tunneling barrier are illustrated: at 
ts the electron is inside the atom in front of the barrier, at 
tr it leaves the barrier at the tunnel exit and subsequently 
moves in the classically allowed region, (b) Real parts of the 
saddle point times for a monochromatic laser with the electric 
field E{t) = Eq cos Lot and vector potential A{t) = —Aq sincot, 
Eq = Aqlj. The superscripts of tr indicate the signs (-|-), ( — ) 
of the asymptotic momentum p^; the numbers 1 and 2 indi- 
cate whether the corresponding quantum trajectory is "long" 
or "short," respectively. The negative vector potential (red) 
and the electric field (blue) in arbitrary units are indicated. 
Panels (c)-(e) illustrate the asymptotes of real trajectories in 
perpendicular direction 6 — tt/2 (c), for an angle < 9 < n /2 
(d), and parallel to the polarization direction, 6 = (e). The 
bold solid lines pg'^ and pg^' represent trajectories in plain 
SFA where initial and asymptotic momentum p are equal. 
Thin solid lines represent the corresponding TCSFA trajecto- 
ries, which start ofi^ like the dashed trajectories with a mod- 
ified initial momentum but then get Coulomb-distorted such 
that the asymptotic momentum is p. 



so that the amplitude ([T]) becomes 

^(TCSFA) ^ ^Cp„(4"'))e^^(*=°''). (8) 



Here, we write Pq to indicate that p ^ const now. Instead 

p{ti" ') ~ Pq. The saddle-point times are also modified 
(hence a'), as they satisfy 

1 
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(9) 



As illustrated in Fig. [TJa), the integration in complex 
time in W(ts) can be splitted into two parts. 



W = W' 



sub 



(10) 



The sub-barrier action W^^^ = — • • • , where tr = 
Rets, may be calculated along the path parallel to the 
imaginary axis down to the real axis. Provided the in- 
tegrand is analytic, the result for the total action W is 
independent of the integration path. The action = 
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— ■ ■ ■ in the classically allowed region is evaluated 
along the real axis. 

The well-established imaginary time method (ITM) 
has been used to study quantum tunneling through 
(time-dependent) barriers in various contexts (e.g., in 
|22|). At t = ts, the electron is in the bound state with 
dynamical variables in general being complex. At t ~ tr, 
the electron emerges at the outer turning point r{tr) 
(tunnel exit). The boundary conditions for the quantum 
orbits may be chosen such that from that time tr on all 
variables are purely real [23j . Let us first summarize the 
Coulomb-correction to this real part of the propagation, 
as it has been employed in 0, [ij, [l^ [l3| . 

A. Real-part Coulomb correction 

For W^'^, the real trajectories with the Coulomb in- 
teraction taken into account are obtained by solving the 
equation of motion ([6]) with the initial conditions 

r{tr) ^ a{tr) -Rca{ts) (11) 
(tunnel exit) where 

a{t) = J A(t')dt' (12) 
is the electron excursion, and 

V{tr) = Po + Mtr). (13) 

The corresponding saddle-point times t'f are deter- 
mined via Eq. ([9]). The equations of motion are solved 
numerically, along with \ up to the end of the 

laser pulse. The asymptotic momentum p'" can then 
be determined via Kepler's formulas. The values of 

•< ' ^ other entities of interest are saved in a data 
base. In the previous works [H, [ij, [l3| the uncorrected 
l^sub,(a ) j-j without Coulomb potential) of the plain 
SFA with p = Pq was evaluated analytically. Finally, 
photoelectron spectra |Mp|^ can be calculated by adding 

coherently all the Mp" ■* of the quantum orbits from the 
data base that have their asymptotic momentum p^" ' 
in a bin around the final momentum of interest p. The 
size of the bin determines the resolution of the photo- 
electron momentum spectrum. In the current work we 
used 1600 X 400 bins to sample the final momentum in- 
terval [—2, 2] X [0, 1] in the spectra of Fig.[5][b)-(d) below. 
Thanks to the azimuthal symmetry with respect to the 
laser polarization direction it is sufficient to consider 
only, e.g., > 0. We covered randomly the initial mo- 
mentum interval [—2, 2] x [0, 2] by 32 million points. For 
each initial momentum several saddle-points contribute. 
In the case of the 3-cycle pulse considered below we thus 
dealt with ~ 100 million trajectories. Because of focusing 
along the polarization axis the randomly sampled initial 
momenta in lateral direction p^o should cover a larger in- 
terval than the final p^ of interest. On the other hand, it 



rarely occurs that \pzo \ > \pz\ so that the initial momen- 
tum interval in polarization direction may be chosen the 
same as the final one. If, for an initial momentum p^o > 
the final momentum p^ turns out to be < we simply re- 
verse the signs of both p^o and p^ . The finer the binning 
and the higher the dynamic range of interest the more 
trajectories are required to ensure good enough statis- 
tics. It is also advisable to shoot randomly in the initial 
momentum plane, not in a regular, grid-like pattern, in 
order to avoid artificial structures in the final momentum 
spectra. The numerical effort to obtain TCSFA spectra 
like the one in Fig. ^c) by solving the equations of mo- 
tion for 100 million trajectories is significant but trivially 
parallelizable (25| . 

B. Sub-barrier Coulomb correction 

One may attempt to solve the equation of motion ([5]) in 
complex space and time under the barrier as well, for in- 
stance, by propagating from the tunnel exit "backwards" 
in imaginary time from tr to ts- However, we found 
that the photoelectron spectra obtained in this manner 
are completely spoiled by artifacts such as caustics, de- 
stroying the interference pattern observed in the ab ini- 
tio TDSE results. If one wanted to propagate in forward 
direction from ts to tr it is not clear what the initial con- 
ditions should be. Choosing, e.g., Rer(ts) = one would 
have to iterate Imr(is) and the imaginary v{ts) such that 
at the tunnel exit position and velocity become real, and 
Pq = v{tr) — A{tr) because tg is determined bypg. In 
the derivation of tunneling ionization rates |18l420| such 
problems are circumvented by matching the action M^'^"^ , 
evaluated with the zeroth order (i.e., plain SFA) quantum 
orbit with the highest statistical weight, to the exponent 
of the unperturbed, initial state wavefunction. This is 
justified, as for r — !■ the Coulomb potential certainly 
dominates the external laser field. Moreover, the quan- 
tum orbits are almost linear (26| and spatially confined in 
the polarization direction, pointing from the origin to the 
tunnel exit. We pursue a similar strategy in this work in 
so far that we consider the plain-SFA quantum orbits for 
the evaluation of W^^^ only. Of course, for the calcula- 
tion of photoelectron spectra we cannot restrict ourselves 
to just the dominating quantum orbits. 

The plain-SFA quantum orbits fulfill 

v{t) = po + Ait) (14) 

where pg is the randomly sampled momentum from the 
real-part propagation described in the previous Subsec- 
tion. The sub-barrier trajectory is 

r{t) = v{t')dt' - v{t')dt', (15) 

which ensures that r{tr) is purely real and r{ts) is purely 
imaginary. We write the sub-barrier action as 

^sub ^ -^sub,0 _|_ i^subcc^ j--|^g-j 
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where 



sub,0 



is the plam-SFA action, and 



-y{t)+ip 



dt 



subc 



z 



■.dt 



(17) 



(18) 



is the sub-CC action. Here, we write on purpose 
Z/^/r^t) instead of the non-analytic Z/\r{t)\. 

In the preceding studies only W = VF^^.o ^ 

W'^° was taken into account whereas in this work, the 
action takes the full form 



W 



sub,0 



subc 



(19) 



including the sub-CC action ^™bcc_ 

The action may diverge if the complex trajectory 
r{t) = 0. However, thanks to the negligible probability 
of sampling momenta pp that lead to such singularities 
they turned out to be harmless in the actual numerical 
implementation of our sub-CC. 



III. INTRA-CYCLE INTERFERENCE 

In order to investigate the influence of the sub-CC 
on the interference pattern in photoelectron momentum 
spectra we first consider H(ls) [Ip = 0.5) in a linearly po- 
larized laser pulse. The vector potential of the laser pulse 
of carrier frequency uj is given (in dipole approximation) 
by 



A{t) = -Aa sin^ ( — ) sinwi 



ujt 



(20) 



for < t < Nc2tt/lj and zero otherwise. Eq = Aquj is 
the peak electric field strength, Nc is the number of laser 
cycles. 

Figure [211 a) shows the photoelectron momentum spec- 
trum, calculated from the exact numerical solution of the 
time-dependent Schrodinger equation (TDSE) using the 
QPROP code [131. The laser parameters are Nc = 4, 
Lj = 0.0228, Eo = 0.0534 (corresponding to A = 2 ^m 
and an intensity of lO^''^ W/cm^). Here we focus on the 
almost vertically aligned interference fringes |28|. These 
fringes stem from the "intra-cycle interference'^!^ [Tsl of 

two trajectories "born" at times ir^'' and tr^^ at the re- 
spective tunnel exits within one laser cycle. The electric 
field has opposite signs at t^^'' and ti^\ meaning that the 
tunnel exits point in opposite directions. In the context 
of the SPA, the two trajectories are typically denoted by 
"short" and "long" according to whether the signs of the 
projections of r (i^) and pg onto the polarization direction 
are equal or opposite, respectively. In the former case the 
photoelectron drifts directly from the tunnel exit to the 



detector whereas in the latter case it starts off into the 
"wrong" direction before it is reversed by the laser field. 

The vertical fringes are also present in the SFA spec- 
trum Fig. [IJb). However, the positions are shifted rel- 
ative to the TDSE result, as indicated by the vertical, 
dashed lines. In fact, when there is constructive interfer- 
ence in the TDSE result, there is destructive interference 
in the SFA spectrum, and vice versa. This indicates that 
there is an erroneous phase-shift between the long tra- 
jectory and the short trajectory of tt in the plain SFA 
as compared to the exact TDSE result. It is expected 
that incorporating the Coulomb-attraction between the 
outgoing electron and the parent ion should cure this dis- 
agreement. 

Let us now "switch-on" the Coulomb-attraction for the 
propagation in the classically allowed region (i.e., from 
the tunnel exit to the detector), as described in Section 
III Al The result is shown in Fig. [5Jc) . As discovered pre- 
viously 0, side-lobes are reproduced in that way. How- 
ever, the vertical interference fringes still have the same 
position as in the plain-SFA result, meaning that the er- 
roneous phase-shift between the long and the short tra- 
jectory of TT is not cured by the Coulomb-correction out 
of the Coulomb-barrier. Only if the sub-CC, as described 
in Section IIIBl is incorporated as well, do the fringe 
positions coincide with the ab initio TDSE result [see 
Fig. [^Jd)]. Hence, the sub-CC not only affects the statis- 
tical weight of a quantum orbit via the imaginary part of 
l^subcc ^^-j- g^jgQ -j-j^g interference pattern in photoelectron 
spectra via the real part of ]4/subcc^ This finding (and its 
analytical derivation in the following) is the main result 
of this paper. 

Let us suppose that only two quantum orbits are dom- 
inating the transition amplitude for a certain asymp- 
totic momentum p. Given the two individual tran- 



sition amplitudes Mp ^ and Mp', the probability to 
find a photoelectron with asymptotic momentum p is 



(2) 



(1) 



(2) 



with M, 



(a') 



e"^^ (a 



1,2). The action I^(° ' is of the form 
P9p . Interference with pronounced contrast requires 
the weights (i.e., absolute values) of Mp^ and Mp'' be- 



ing comparable, e 



The assumption 



that this is the case then yields an ionization proba- 
bility w(p) ~ 2e"^^™^'"(l + COS0). Hence, the inter- 
ference pattern is determined by the phase difference 
(j) = Re[M^(^' — For cj) = mn with m even we 

have interference maxima, for m odd interference min- 
ima, provided the prefactors C^(ti"'' ) in ([1]) do not intro- 
duce extra phase-differences |29| . 

The separation of the action (|19p suggests that we split 
the phase <f> accordingly, 



iSub.O 



/ subcc 



(21) 



From the simulation results in Fig. [5] we know already 
that it is 0*5"*''='= that is responsible for the intra-cycle 
fringe shifts towards the correct positions, as predicted by 
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Figure 2. (color online). Logarithmically scaled photo- 
electron momentum distribution in the pz-px plane {pz = p\\ 
along the polarization direction) calculated using (a) TDSE, 
(b) plain SFA, (c) TCSFA without sub-CC, and (d) TCSFA 
with sub-CC. The calculation is for H(ls) in a 4-cycle, sin^- 
envelope, near-infrared laser pulse of wavelength A = 2^m 
and intensity 10^* W/cm^. Four orders of magnitude in prob- 
ability are shown. 



the ah initio TDSE result. In the following we will derive 
an analytical expression for (jf^^'^'^ in the long- wavelength 
limit. 



A. Fringe shift in polarization direction 

In the long-wavelength limit (oj 0) it is possible to 
evaluate the sub-CC phase difference 0'^"^)'='= in the polar- 
ization direction analytically. Consider a monochromatic 
laser field 



E{t) = Eq coswt, A{t) ^ - A 



sinwi 



(22) 



where Eq — Aquj, and we assume Aq > 0. We will prove 
in this Subsection that in the long-wavelength limit we 
obtain a momentum-independent sub-CC phase shift of 



10^ 



(23) 



In particular, for the case of Fig. [2] where Z = 2Ip = 1 we 



have |(/)' 



subcc 



turning constructive interference max- 



ima of the plain SFA into destructive interference minima 
in the TCSFA with sub-CC and vice versa. As we know 
already from the results in Fig. [2] that the Coulomb cor- 
rection of the propagation in real space is not responsible 



for the fringe shift, we ignore its effect on Pq^-* and 



and thus also on ti^' and ti^^) in this Section and assume 



(aii ( 

m 



(2) 

■Po 



P- 



(24) 



For a monochromatic laser field (P^ . p = {px,Pz), and 
tg — t,- + iti the saddle point equation ([5|) yields 



sin (wtr) cosh (wti) — 



(25) 



cos [iotr) sinh (wtj) 



T^'^V^. (26) 

^0 



Equation (|25p shows that the possible tr can be deter- 
mined as the intersections between p^J cosh(a;ti) 2± pz 
and the negative vector potential —A{t) = Aosinwt, as 
indicated in Fig. [IJb) . The figure shows that a pair of 
saddle points per cycle, e.g., t^^^^ and t^^~^\ exists for 
a specified Pz (dashed horizontal lines). The so-called 
intra-cycle interference is thus due to two saddle-points 
for which 



i(2) 



P'=P'=P, 
fW - _f(2) 

.(1) 



UJ 
(2) 



(27) 



t. 



t, 



are fulfilled. In the polarization direction, px = 0, the 
sub-CC action is 



subc 



■.dt = 



dr (28) 



with z{t) = z{tr + it) along the integration path t = 
tr -I- ir, T G [ti , 0] , as shown in Fig. [TJa) . The superscript 
(a) is suppressed here. For the sub-barrier trajectory in 
polarization direction calculated from (jl5[) 



z{t) 



An 



[coswi — Re coscl><s] + Pz[t — Rctg] (29) 



results. Separating the real and imaginary parts of z{t), 
z{t) = a(r) + ih{T), (30) 

one finds 

a(T) = — -cosujtj. (coshwr — coshwfi) , 

A, (^^^ 
6(t) = PzT ainujtrSmhujT. 

UJ 

In the long wavelength limit w — > 0, the functions sinh 
and cosh in a(T) and 6(t) can be expanded about w = 0, 
leading to 



a(r) 
6(r) 



Ao 
2 

tA 



U) cos (jjtr 



(32) 



with 



A ~ — Aq sin ujtr 



—^ui^tfsva.ujtr- (33) 
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Let us consider emission times tr for positive p^. In is relevant for interference patterns in the spectra. In- 
this case < ojtr < tt and thus A > 0. As t > serting ([5^ leads with for = and w — >■ 0, 
[3l| as well, 6(t) > 0. Hence the integrand in 
p^subcc (j^gj) jnay be simplified [by choosing the pos- 
itive root] to Z [P{t)]~^^'^ = Z/S{t) for a(T) > and 

Z [z^{t)] ^^'"^ ^ -Z/S{t) for air) < so that it is suffi- 
cient to consider the simplified action 



~z{r) 



The real part of this action, 

rSubcc' 



ReWr-^±l ^^.(^^^r, (35) to 



Aouiti cos ujtr ~ =F-\/2/p, (36) 



ReW: 



subcc 



±- 



2Zt^A 



■dT 



= ±- 



arctan 



^2lp - A2 At2 




(37) 



It follows from ((5^ and ((77|) that, for a given p^, the two 

saddle-point times t'^p and ii^'' are such that a^^^r) = 
—a^^^r) and 6^^'(t) = 6(2)(t). The opposite signs of 
a(T) have been considered already in our calculation of 



jsubc 




(38) 



as claimed in ([221) • 

For a hydrogen-like ion in a state with principal quan- 
tum number n = 1, 2, 3, . . . we have Ip = Z'^ /(2n'^) and 
thus 1 1/)^^'"^'^ I = mr. In particular, for ground states 7i = 1 
we have independently of Z that 1 1/)**"'^'^'^ | = tt. This is 
the required phase shift for bringing the intra-cycle fringe 
pattern in agreement with the ab initio TDSE result [see 

Fig.[ird)[. 

We compared TDSE and TCSFA results for other pa- 
rameters Z, Ip, and n and find good agreement if the sub- 
CC is included, provided the minimum number of pho- 
tons required for ionization Ip/uj 3> 1 and the Keldysh 
parameter -y/Jj,/ {2Up) ^ 1, where Up = Aq/ 4 is the pon- 
deromotive energy. Moreover, the smaller Z the better is 
the approximation to use the zeroth-order, plain-SFA tra- 
jectory for the sub-CC. Excited states n > 1 may intro- 
duce extra phase-differences via the prefactors Cp{t'f ^) 
in ([1]) |29]], not discussed in this work. 

The reader may have noticed that the precursors of 
ATI rings, visible in the exact TDSE resuh of Fig. [D(a) 
and in the SEA result (b) as circular structures around 
the origin, are actually less clear in the TCSFA (c), let 
alone in the sub-CC TCSFA (d). The reasons for this 
deterioration are the overemphasized classical boundaries 
and caustics that spoil for few-cycle pulses the low-energy 
part of the spectra by rendering the few trajectories that 



contribute to the ATI-ring structure relatively less im- 
portant. The situation improves for more-cycle pulses, 
as in Fig. [3] below, where more trajectories contribute to 
inter-cycle features such as ATI rings. 



B. Interference in perpendicular direction 

We are also interested in effects from the sub-CC in 
the perpendicular direction pz = 0. In the SEA, the 
short and long trajectories are asymptotically parallel, 
as shown in Figs. [IJc-e). In the perpendicular direction 
[Fig.[IJc)[, the parallel asymptotes are vertical when two 
saddle points merge as in ([^5)1 . One easily ver- 

ifies that in Fig. \V[h) ti^^^ merges with ti^ ■* and t^'^'^ 

merges with ti^ \ Hence, in plain SEA the momenta 
and the saddle point times of long and short trajectories 
become degenerate in perpendicular direction [Po^^ and 

Pq in Fig.[llc)[. However, with the long-range Coulomb 
interaction included, this degeneracy is lifted. In order 
to have a common asymptotic momentum p' the initial 
momenta pg must be different, as indicated in Fig. [IJc). 

In order to cancel the longitudinal shift induced by 
the long-range Coulomb distortion, the two trajectories 
should satisfy 



POz 



(2) 
POz ' 



POx 



(2) 
POx- 



(39) 



The corresponding saddle points are easily found for 
these criteria from Fig. [T](b). Taking, e.g., the saddle 
point (1+) as the reference, the other saddle point is 
(1 — ). Therefore, the initial conditions for these trajecto- 
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(1+) 



Pqz 



(1+) 



POx I 



4i+) = ^ + t(i- 



(40) 



In the following, the phase contributions 4>l"^'^^', (^^']^'° 
and are deduced, and the interference pattern in 
the perpendicular direction is inferred. The subscript 
"c" indicates that the phase difference between the tra- 
jectories has been generalized to allow for the real-space 
trajectory-distortion by the Coulomb force. The interfer- 
ence pattern in perpendicular direction is determined by 
the symmetry properties of the interfering trajectories, 
while it is independent of the detailed. Coulomb- modified 
dynamics because both trajectories of a pair collect equal 
phase distortions. This suggests that the SFA, TCSFA 
with or without sub-CC result in the same interference 
structure in the perpendicular direction, as will be dis- 
cussed in the following. 



1. Phase difference from sub-CC: 



i subcc 

Pc± 



We show first that the sub-CC phase difference (f>l^'^ 
does not contribute in (|2T]). 



0; 



c± 



Re 



W 



subcc(l + ) 



Libcc(l — ) 



= 0. 



(41) 



In passing, we would like to mention that caution has 
to be exercised regarding which saddle-points are chosen 
for = if zeroth-order, plain SFA were used, for which 
the degeneracy of the saddle points is not lifted. If the 
"usual" long and short trajectories (1+) and (2+) with 

ti^^^ and ir^^'' shown in Fig. ]V[h) are chosen, z^^+^ (r) 
can be expressed by variables of (1+) as 



5(2+) 



Ao 



cos{ujtl,^^^) (coshwT — coshujti) 
i^'^'^T sinwtl^'*''' sinhwT 



(45) 



Compared to (|43p . this expression differs by a sign only 
in the first term. Since i^"'^"''-' (r) = i'-^+^(T), the imagi- 
nary parts of z{t)'^ +x{t)'^ in the denominator of (|^^ for 
(1-I-) and (2-|-) are opposite in signs. It seems a paradox 
occurs: in the exact perpendicular direction, integrands 
are identical since the imaginary part in (j45|) is zero, and 
thus no phase difference exists. But a tiny longitudinal 
momentum pz would lead to the abrupt appearance of a 
finite phase difference, which is unphysical. The problem 
does not arise in the above employed sub-CC, thanks to 
the Coulomb distortion of the trajectories in real-space, 
which enforces Pq^^'' = — Pqz ^ thus determining the rel- 
evant saddle-point times to be (1-I-) and (1 — ), not (1+) 
and (2+). 



The real part of the sub-CC action reads 

dt 

It^ ^z{t)^ + x{ty 



ReW 



subcc 







dT 



U ^f(T)2 +i(T)2 



(42) 



For a monochromatic laser pulse (P^ . the zeroth-order 
sub-barrier trajectory (jl5|) as a function of t = tr + ir 
can be cast into [see (l30l). (1311)1 



5(t) 



An 



COS (iutr) (coshwr — coshcjt,;) 



+ « PzT 



A„ 



sinwtr sinhwr 



(43) 



and 



xir) = ip^T. 



(44) 



Substituting tr^^-* and pi^^"* into z{t) (|43|) . we find with 
the help of the conditions ^ 5^1+) (r) = -z'^^-\t). 
For x{t) ([33]) clearly follows x(^+)(t) = ^(^"^(t). As a 
consequence, the integrand in ([31]) is identical for both 
trajectories. As the integration starts for both trajec- 



,(1+) 



tories at t 
and (|4ip is proved. 



(1-) 



ti we obtain Re W[ 



subcc( 1 + ) 



Re w^"'^'='=(i )_ i ^}^g sub-CC phase difference vanishes. 



2. Phase difference from sub-barrier propagation without 
Coulomb interaction: 61^}'''^ 



After having identified the relevant saddle-points (1-f ) 
and (1 — ) we now show that 



,sub,0 



= Re 



W 



sub, 0(1+) 



w 



sub,0(l-) 



= 0. 



(46) 



Substituting a monochromatic laser field (P^ into (flT)) 
yields 



^sub.O _ ■ I r , 



= i{Ip + -p' + Up) ti 



AqPz 

ijj 
Al 



(coswir — cos wis) 



(47) 



H 2. (sin2ajtr ^ sin2ajts) 

where p = \/p\ + p^. The real part reads 



ReP^sub,o 



coswir(l — coshwti) 



AqPz 

— - sin 2u}tr(l — cosh2a;ti). 
8cj 



(48) 



Inserting pi^^'^ and t^^^\ which satisfy (PO)) . we find 



j^^^sub,0(l-) ^ j^^^sub,0(H-) 



(49) 



confirming (|46p . 
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3. Phase difference from real propagation: (pl"^ 

The interference pattern in the perpendicular direction 
is entirely determined by the real-time trajectories and 
the corresponding phase difference 

Given real trajectories z{t), x{t), Vz{t), Vx{t) for t £ 
[triTp], the action due to the propagation in real time 
is 



dt. 



The symmetry properties for two interfering trajectories 
in the perpendicular direction are 



,(1 



(t)=-.r)(t+^ 



(50) 



^^"^ ^ -tt/w. Casting Wl 



Moreover, tr 
of the variables of (1+) one finds 

,(1+) \2 
Z 



in terms 



W 



rc(l- 



Tp+Tr/cj 



.(1 + ) 



The phase difference (ff^^^ 
given by 



W 



rc(l + ) 



J- 



M^r^'-) is thus 




I{t) dt= I I{t) dt, 

J Tp 

(52) 

independent of the specific saddle-point pair under con- 
sideration. As Tp — > oo the asymptotic position in po- 
larization direction remains finite, 2(t)|t_i.oo ~^ Zoo, but 
x{t)\t-foo ^ oOi leading to a vanishing phase difference 
due to the Coulomb potential. In the kinetic energy term 
of the integrand, ^^(Olt^oo M^) and Vj:{t)\t^oc Px, 
so that 



-pl + -A^{t) + Ip ] dt 



1 



(53) 



J^(t) dt. 



Using -px = ■\/2 (nuj — Up — Ip) from the energy 
conservation for ATI in the perpendicular direction 
with n the number of absorbed photons and Up = 

^2 27/:; lT^^A\t)dt we and 

€1 = nn. (54) 

This result agrees with the findings in Ref. Q that every 
other ATI ring is destructively "interfered away" in per- 
pendicular direction. Because of the pre- factor Cp{ti°'^) 
in ([1]) the interference is destructive for odd n in the case 
of even-parity initial states and vice versa. 

The interference concerning the ATI rings in perpen- 
dicular direction predicted by ([M)) is confirmed by the 
results for ionization of II(ls) in a 9-cycle trapezoidal- 
envelope laser pulse of 800 nm wavelength shown in 
Fig. [31 The other laser parameters are given in the figure 
caption. 



IV. ATI ON-RING STRUCTURE 



ATI peaks appear as rings of radius 



P = y^2 {nuj - Up - Ip) 



(55) 



in the photoelectron momentum spectra. It is known 
from experiments that the ionization probability along 
the ATI rings is not isotropic but exhibits nodal struc- 
tures. These structures, also visible in Fig. [31 have been 
analyzed recently within the SFA Q. The radius of the 
ATI rings is well reproduced already by the plain SFA, 
which is verified by comparing Fig. with the ab ini- 
tio TDSE result Fig.[3I^a). However, the positions of the 
nodes along low-order ATI rings do not agree well. Not 
even the number of minima and maxima agrees. Take, 
e.g., the second ATI ring {p = 0.3939, absorption of 14 
photons) . The TDSE result in (a) shows 9 maxima along 
that ring, the plain SFA result in (b) only 7. The real- 
time Coulomb-corrected SFA in (c) does not affect the 
number of maxima. Instead, the sub-CC SFA in (d) does 
give the correct number. Along the next ATI-ring (ab- 
sorption of 15 photons) TDSE and plain SFA agree al- 
ready (10 maxima), and TCSFA with or without sub-CC 
does not change the nodal structure. Along the first ATI 
ring the TDSE shows 8 maxima, which is not properly 
reproduced, even with the sub-CC TCSFA, probably be- 
cause our sub-barrier Coulomb correction involving only 
the zeroth-order plain-SFA trajectories is insufficient. 

Figure [31 shows the momentum distribution along the 
second ATI ring for the ab initio TDSE, the SFA, and 
the TCSFA with or without sub-CC. TCSFA with sub- 
CC agrees well with the ab initio TDSE result. The 
results from the TCSFA without sub-CC and the plain 
SFA are similar. They agree well with the TDSE around 
9 = 71/2 but deviate for smaller angles. As discussed 
above, there is one maximum less in the plain SFA and 
the TCSFA without sub-CC, as compared to the TDSE 
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-2.0 -1.5 -1.0 -0.5 D.O 0.5 1.0 1.5 Z.O 



Figure 3. (color online). Same as in Fig. [2] but for a 9-cycle, 
trapezoidal-envelope laser (2, 5, and 2 cycles for up-ramping, 
constant-envelope, and down-ramping, respectively). The 
wavelength is A = 800 nm, the intensity is 10^'' W/cm^. In 
the spectrum calculated using the TCSFA with sub-CC (d) 
only trajectories of type Tl and T2 in the nomenclature of 
fv^ have been considered (3^ . 




Figure 4. (color online). The ionization probability along the 
second ATI ring with p = 0.3939 in Fig. [3] according to the 
ab initio TDSE, the (plain) SFA, the TCSFA (without sub- 
CC), and the TCSFA with sub-cc result. The polar angle 9 
is in the range [0, 7r/2], 6 = corresponding to px = 0, and 
6 — ir/2 corresponding to p, = 0. The plain SFA result is 
as noisy as the TCSFA because it was calculated using the 
same statistical sampling approach as for the TCSFA (which 
is explained briefly in Sec. Ill Al and in more detail in ]V^). 



V. SUMMARY 

We extended the previously introduced trajectory- 
based Coulomb-corrected strong field approximation to- 
wards the inclusion of sub-barrier Coulomb effects. To 
that end we used the plain strong field approximation 
quantum orbits and inserted them into the Coulomb- 
corrected sub-barrier action. The (additional) phase- 
difference between the so-called "long" and "short" trajec- 
tories introduced in that way has a pronounced qualita- 
tive effect on the photoelectron spectra, namely a phase 
shift of Ztt/ along the polarization axis for linear 
polarization in the long- wavelength limit. In the case 
of ionization from the ground state of atomic hydro- 
gen or hydrogen-like ions, the phase shift is tt. turn- 
ing constructive intra-cycle interference maxima into 
destructive interference minima and vice versa. We 
benchmarked the predictions of the trajectory-based 
Coulomb-corrected strong field approximation including 
sub-barrier effects by comparisons with ab initio solu- 
tions of the time-dependent Schrodinger equation. We 
also showed that the sub-barrier Coulomb correction af- 
fects the nodal structure along low-order above-threshold 
ionization rings. 

Our semi-classical method is applicable to sufficiently 
long wavelengths (/p/w ^ 1). However, the nuclear 
charge Z (and thus Ip) has to be small enough to render 
the zeroth-order orbit a good approximation for the sub- 
barrier Coulomb correction. Moreover, the trajectory- 
based Coulomb-corrected strong field approximation dis- 
plays too many unphysical caustic structures in the spec- 
tra if the Keldysh parameter \/Ip/ (2[/p) is not small 
enough (i.e., around 1 or smaller). It thus turns out 
that our method is applicable in a regime where the di- 
rect ab initio solution of the time-dependent Schrodinger 
equation in position space is stretched to its limits: long 
wavelengths and high ponderomotive energy. 

Our method is, in principle, also applicable to more 
complex systems for which strong field approxima- 
tions already have been or can be developed. The 
most straightforward extension is towards two Coulomb- 
centers (e.g., 112"). Also non-sequential ionization in ef- 
fective two-electron systems is accessible. However, due 
to the higher-dimensional momentum space to be sam- 
pled in such systems it will be numerically very demand- 
ing to obtain differential photoelectron spectra with good 
enough statistics. The infiuence of many-electron polar- 
ization effects could be investigated by incorporating a 
polarization potential into the equations of motion for 
the quantum trajectories (sslls^. 



and the TCSFA with sub-CC. Hence we conclude that 
the nodal structure of low-order ATI rings is affected by 
the sub-CC. Higher-order ATI rings may be effected by 
the side-lobes (clearly visible in Fig. ISUa) and discussed 
mi). 
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